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Parallaxes of W UMa stars in the Hipparcos catalogue have been analyzed. 3 1 W UMa stars, which have the most accurate 
parallaxes {(JT,/Ti < 0.15) which are neither associated with a photometric tertiary nor with evidence of a visual compan- 
ion, were selected for re-calibrating the Period-Luminosity-Color (PLC) relation of W UMa stars. Using the Lutz-Kelker 
(LK) bias corrected (most probable) parallaxes, periods (0.26 < P{day) < 0.87), and colors (0.04 < {B ~ V)o < 1.28) 
of the 31 selected W UMa, the PLC relation have been revised and re-calibrated. The difference between the old (re- 
vised but not bias corrected) and the new (LK bias corrected) relations are almost negligible in predicting the distances 
of W UMa stars up to about 100 parsecs. But, it increases and may become intolerable as distances of stars increase. 
Additionally, using (J — H)o and {H — Ks)o colors from 2MASS (Two Micron All Sky Survey) data, a PLC relation 
working with infrared data was derived. It can be used with infrared colors in the range —0.01 < {J — H)o < 0.58, and 
—0.10 < {H — Ks)o < 0.18. Despite 2MASS data are single epoch observations, which are not guaranteed at maximum 
brightness of the W UMa stars, the established relation has been found surprisingly consistent and reliable in predicting 
LK corrected distances of W UMa stars. 

© 2008 WILEY- VCH Vcrlag GmbH & Co. KGaA, Wcinhcim 



00 

a\ 
^. 

i> 

o 
oo 

o 



X 



1 Introduction 

Low-mass contact binaries, popularly known as W Ursa Ma- 
joris (W UMa) stars are easy to recognize by a light curve 
with equal (or nearly equal) depths of minima which are 
wide enough to touch one another. Even if the total lumi- 
nosity produced solely by the more massive component, 
efficient energy distribution through a common envelope 
makes the their sur face brightness practic ally the same over 
the visible surface (lRucinskilll985ill993h . Even though the 
effective temperature is the same on the surfaces of com- 
ponents, there could be different masses hidden under the 
common atmospheres. The ir mass-ratios span from almost 



of geometry, Kepler's third law and main-sequence rela- 
tionships. The calibration, which is in the form of My = 
A/y(logP, -B — V), has been widened later to include a 
metallicity term [Fe/H] ( Rucinskill 19951) . Those earlier cal- 
ibrations were based on 18 systems, mostly members of 
high galactic latitude open clusters and of visual binaries, 
including only three nearby W UMa systems with known 
trigonometric parallaxes, 44i Boo B, VW Cep and e CrA. 

After Hipparcos parallaxes became available, the Period- 
Luminosity-Color (PLC) relation has been re-calibrated with 
40 nearby W UMa systems in the solar neighborhood with 
parallaxes having relative errors from 2.7 to 24%. With this 

new data, the accuracy of the calibration (±0.25 mag at 

unity (V753 Mon, q=0.973. lRucinski. Lu & Mochnacki','2000') a predicted My) was claimed: "it is no longer limited by 
to ver y small values as small as g=0.066 (SX Crv, Rucinski et althe parallax data but, paradoxically, by the lack of reliable 
200lb . Therefore, W UMa binaries are non-equilibrium sys- photo metric da ta" (Rucins ki & Duerbeckl [ l997h . However, 



tems exchanging mass and energy between the components. 
Hence, despite such an external simplicity, their internal 
structure is rather complex. 

External simplicity and the observational properties of 
W UMa systems were the basic re asons for establi shing 
an absolute magnitude calibration by iRucinskil (1 1 9941) . The 
calibration uses two observational quantities, the orbital pe- 
riod and intrinsic color, e.g. {B — V) or {V — I). The pe- 
riod and the color are correlated through a combined effect 
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later, iRucinskil (l2000t i2004l) argued: at earlier times, when 
limited information w as available, it se emed that metallicity 
dependence did exist dRucinskil Il995h . But later when ex- 
tensive data for many globular clusters we re examined, th e 
metallicity term [Fe/H] become obsolote (IRucinskil 120001) . 
At least, solar neighborhood calibration works well for con- 
tact systems of different metallicities at level of uncertainty 
(7 = ±0.28 mag. Thus, ther e is no need to k eep the metal- 
licity term in PLC relation (iRucinskil |2004|) . At last, after 
metallicity term was dropped, the color term was also re- 
moved, and the re l atio n is transformed to a much simpler 
form by JRucinskii (|2006|) . when estimating spatial density 
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of contact binaries based on the ASAS survey as a necessity 
because the lack of color indices in the ASAS data. This cali- 
bration was established from 21 systems of good Hipparcos 
{a-n/n < 0.12) data which were carefully chosen as being 
free of triple and multiple systems. 

PLC relation provides absolute magnitudes of contact 
binaries if their de-reddened colors and periods are known. 
Consequently, PLC relation has many practical usages. It 
is useful not only in calculating space velocities and space 
densi t ies of W UMa (JBilir et aliboost lOettel. Geske & Mckav , 
2006t iRucinskil 120021) stars but also for co nfirming or dis- 
provi ng of contact binaries in stellar clusters ( IRucinskil Il998 , 
20001) . Moreover, it permits W UMa stars to be used as stan- 
dard candles. In fact, through PLC relation, W UMa stars 
become more likely to be used as standard candles than RR 
Lyr stars because W U Ma stars are 24000 times mor e com- 
mon than RR Lyr stars Wucinski & Duerbeckll 1997b . 

Unfortunately, random errors, presumably symmetric on 
the measured parallaxes, do not provide symmetric uncer- 
tainties on the computed distances. Therefore, a measured 
trigonometric parallax is very likely to be larger than the 
true paral lax. The problem has already been noticed and 
studied by Lutz & Kelken(ll973b . Assuming a uniform space 
distribution of stars and a Gaussian distribution of observed 



parallaxes over a true parallax, iLutz & KeUcen (Il973h have 
revealed that there exists a systematic error in the computed 
distances which depends only upon the ratio cttt/tt, where 
TT is the observed parallax. However, this crucial study has 
not produced proper response in the literature. Many studies 
including existing PLC calibration we re completed withou t 
even discussing the Lutz-Kelker bias. Ijerzvkiewicj (120011) 
has shown that only the careful studies which used paral- 
laxes with ct^/tt < 0.1 could be excused since the bias 
would be negligible. Otherwise, not taking the Lutz-Kelker 
bias into account would either alter the conclusions or inval- 
idate them altogether if relative larger errors were involved. 

Standard L utz-Kelker corrections b ecome significant if 
ct^/tt > 0.05 (Maiz Apellaniz, '2005). Therefore, Period- 
Luminosity relation of .RucinskL (,2006,) . which relies on 21 
W UMa stars with sufficiently accurate {ut^j-k < 0.12) par- 
allaxes, may not be excuse d. The classical PLC relation of 
Rucinski & Duerbeckl (11997.) . which was founded on paral- 
laxes of 40 W UMa stars with relative errors from 2.7 to 
24%, definitely needs to be re-calibrated. This study aims 
to improve existing PLC relation by not only correcting it 
according to Lutz-Kelker bias but also refining the sample 
of W UMa stars by eliminating the ones which have less re- 
liable parallaxes because of belonging to multiple systems. 

Moreover, applicability of the PLC relation at infrared 
colors are investigated. Being based on similar principles, 
a similar relation (PLC) for cataclysmic variables (CV) has 
already been established by using (J — H) and {H — Kg) 
colors from 2MASS (Two Micron All Sky Survey) data and 
proved to be useful for estimating CV distances bv lAk et al.l 
(l2007h . Despite, 2MASS data are compiled from single epoch 
observations, which are not guaranteed to be at maxima, the 
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Fig. 1 Relative error {aT^/ir) histogram of Hipparcos par- 
allaxes for W UMa stars where a) shows cumulative sample 
and b) relative numbers. Vertical lines mark 50% and 15% 
relative errors. 



relation which uses [J—H) and {H — Kg) colors of 2MASS 
photometric system have been found to be as useful as the 
existing relation which involve {B — V) and {V — /) colors 
(IRucinskil |2004|) . The advantage of this new relation is that 
it is less effected by de-reddening problems. 



2 Data 

There are 75 1 W UMa binaries listed in the revised edition 
oftheGCVS. Only 144of them were found to have trigono- 
metric parallaxes in the Hipparcos catalogue (iESALIl997l) . 
Nine of them could be discarded right away because their 
relative errors (ct^/tt > 1) are intolerable. For the rest (135 
systems), the distribution of relative parallax errors have 
been shown in Fig. 1. 

Calibration process of PLC relation requires both accu- 
racy of calibration data and statistical significance. While 
statistical significance requires as many as possible W UMa 
stars, accuracy of input data, on the other hand, demands 
using only the data of W UMa stars with the most accu- 
rate parallaxes. Therefore, an optimized upper limit for the 
relative error of parallaxes needs to be determined. 

2.1 Lutz-Kelker bias corrected parallaxes 

It is known that a standard Lutz-Kelker (LK) correction < 
AMv >, which corrects mean absolute magnitude of a 
sample, is significant if (j^/tt > 0.05 and increases as a-^/Tr 
increases. Corrections were cla imed to diverge and became 
incomputable if ct^/tt > 0.1 75 dBrow n et al.Ul997!: 



1973T[l9 75; Maiz Apellan i3.l2005llN orgaard-Nielse] 



Lutz & Kelkei . 



Oudmaiier. Groenewegen & Schriiveilll998l : ISmit] 

2003h . Standard LK bias primarily computes an expected 



senljl978l: 
:hr il987alli 



a 



value of parallaxes, which stands for the mean of the true 
parallax distribution of an uniformly distributed sample of 
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stars. Here, we are interested in the most probable value of 
the true parallax (ttq) for a single star which has an observed 
parallax (tt) and associated error (aTr)- As for estimating dis- 
tances, one can transform the posterior probability of the 
true parallax to that of the distance and use it to find the 
most probable distance ( ISmithlll987al) . Considering that the 
present W UMa sample is well within the Galactic disc (TV 
Pic has the biggest z = 127 pc, where z is the distance from 
the galactic plane), it can be assumed that within the Umited 
space at the vicinity of the star, th e distribution is uniform. 
Consequently, according to lSmithI (Il987a) 



Table 1 Selection of the calibration sample from the list 
of W UMa stars with parallax errors less than 17.5%. 



1 



na = ni- + -^T^Wa^J^ 



(1) 



is the relation between the most probable true parallax (com- 
puted) and the observed parallax (tt) and it's associated stan- 
dard error (o-tt)- 

Although Eq. 1 is valid mathematically for a-^/TT < 
0.25, only the most probable true parallaxes of the systems 
with cttt/tt < 0.175 (standard LK limit), and correspond- 
ing distances (do = I/ttq) were computed and displayed in 
Table 1. 

2.2 Other observational data 

In addition to LK bias corrected parallaxes, the most reli- 
able orbital periods and intrinsic colors are needed to re- 
vise the PLC relation. Orbital periods and spectral types 
are mo stly from the catalog of iPribuUa. Kreiner & Tremko 
(l2003h . while the maximum visual brightnesses (Vmax) were 
taken from[Pribulla & Rucinski (2006). Intrinsic {B - V)o 
colors of main sequence stars were taken from Neill Reid's 
WEB page [J. In addition to these visual photometric data, 
the infrared brightnesses J, H and Kg magnit udes were 
taken from the Point-Source Catalogue and Atlas JCutri et al 



2003l:ISkrutskie et al.L 12006 ) which is based on the Two Mi- 



cron All Sky Survey (2MASS) observations. The 2MASS 
photometric system comprises Johnsons J (1.25 /im) and 
H (1.65 /im) bands with the addition of Kg (2.17 /im) band, 
which is bluer than Johnsons X-band. Infrared data are for 
to establishing PLC relation at infrared wavelengths. 

The color excess E{B — V), although a relatively small 
quantity, it can effect the calibration in a systematic way. It 
has already been noticed by iRucinski & Duerbeckl (119971) 
that the col or excess were hkely overestim ated in previ - 
ous studies dRucinski & Kaluznvl 1 1981; Ru cinskil Il983h . 
Therefore, we have carefully re-investigated the color ex- 
cess of the stars in Table 1 using three independent meth- 
ods. The first method uses spectral types to estimate intrin- 
sic color of a system. The second method uses the color 
excess found from the literature directly, which are the esti- 
mates from the hydrogen column density or main-sequence 
fitting. The th i rd me thod computes the color excess from 
Schlegel et al.l (119981) maps by using NASA Extragalactic 
Databascl. Since sample W UMa stars are relatively nearby. 
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a multiple system. 



the color excesses according to ISchlegel et al.l (119981) need 
to be reduced. First, Eoo {B — V) color excess in the galactic 



latitude (6) and longitude (Z) was taken from Schlegel et al. 
([1998). Then, the total absorption towards the star in the 
galactic disk in the F-band was evaluated as 



A^{h)=i.lE^{B-V), 



(2) 



where the subscript symbolizes up to infinity but actually it 
is up to the edge of our galaxy in the line of sight. Thirdly, 



www.an-journal.org 



© 2008 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinheim 



Z. Eker et al.: New absolute magnitude calibrations for W Ursa Majoris type binaries 



0.2 



0.15 - 



0.1 



0.05 - 





' ' 




' ' 


' 


A 




' 






- 


















o 




- 






. * 








X 






X 


- 














O 






o 








oo 


o 


o 
o 


gD 





o 


X 




t 








o 




a 


o 








- 




8 

t 


o + 




X 


xx 




+ 




X 



50 



100 
Distance (pc) 



150 



200 



Fig. 2 Color excesses; (o) observed minus intrinsic, (+) 
literature and (x) Schlegel et al. (1998) reduction. 



the interstellar ab sorption up to the star dist ance d is calcu- 
lated according to lBahcall & Soneiral (Il980l) 



2.3 Formation of the final calibration list 



Regarding the accuracy, today's technology is capable of 
determining a visual brightness within the accuracy of few 
millimagnitudes. Because of O'ConnelFs effect, which could 
be several millimagnitudes dDavidge & Milonel Il984l) . the 
maximum visual brightness (Vmax) of W UMa stars have 
additional uncertainties. Error contribution of a visual bright- 
ness, even if it is in the order of several millimagnitudes, is 
negligible compared to the other contributing errors. For ex- 
ample, uncertainties of E{B — V), which contributes to the 
uncertainty of interstellar absorption th roughly = 3.1E{B— 
V), are in the order of 0.01 -0.02 mag (Rucins ki & Duerbeck , 



A4b)^A^{b) 



1 



exp\ 



d sinih) 



(3) 



where H is the scaleheight for the in terstellar dust which is 
adopt ed to be lOOpc as usual (see e.g. lMendez & van AltenaL 
19981) . Finally, the color excess for the star at the distance d 
is estimated as 



Ea{B-V) = Aa{b) /i.l. 



(4) 



It is demonstrated (Fig. 2) that the first method usually over- 
estimates the color excess with respect to the other two meth- 
ods. The color excesses from the literature and the values 



computed using lSchlegel et al.l (11998 ?) extinction maps agree 
with each other. Literature values could not be trusted espe- 
cially if they are zero, which could be the cases because the 
color excess was ignored b lindly. On the other ha nd, values 
computed by the method of lSchlegel et al.l (Il998h extinction 
maps are not trusted at low galactic latitudes (|6| < 10°) due 
to local inhomogeneities of galactic disc. Therefore, care 
was given to E{B — V) values of all stars one by one and 
the most trustable value was selected. Adopted values are 
usually the ones which are taken from the literature where 
the source references are given or the method of inquiry is 
indicated (Table 2). 

Once color excesses are known, then the intrinsic colors 
and the interstellar absorption in the y-band were computed 
as 



{B ~V)a = [B -V) - E{B -V), 
Av ^i.lE{B^V). 



(5) 
(6) 



Consequently, de-reddening of t he in frared bands becomes 



ready usi ng the relations given bvlFior ucci & Munml ( 2003h. 
(see a lso, lAk et al.Ll2007l:lBilir. Guver & Aslanll2006l:lBilir et 



1997; Rucinski, 2002, 2006). The uncertainty of Ay, how- 
ever, is not only due to the uncertainty ofE{B — V), but also 
due to the uncertainty of the coefficient in front of E{B — 
V). Depending upon the position of a star relative to the 
Sun in the galactic plane, there could be cases where the 
coefficient could be as small as 2 .75 and as big as 3.52 
JCardelli. Clayton & Mathisl 119891) . Therefore a rough es- 
timate of the uncertainty of Ay is in average 0.06 magni- 
tudes. Even if the error contribution of Ay is overestimated 
to be 0. 1 mag, it would still be smaller than the uncertainty 
contribution of a 10% parallax error, which is 0.217 mag. 

There are 29 W UMa systems with ct^/tt < 0.10 in Ta- 
ble 1 , which could be con sidered statistically s i gnific ant at 
least at the limit. However, Prib uUa & Rucinskil ( 120061) have 
found that up to 59% of W UMa stars have companions. Ac- 
tually, one of the criteria of multiplicity for lPribuUa & Rucinski 
(2006) was the large error of the observed parallax. Trigono- 
metric parallaxes are frequently wrong for these systems be- 
cause complexity of the multiple orbital motions shows it- 
self both at proper motion and parallax measurements. There- 
fore, W UMa binaries which are known to be a ssociated 
with companions were excluded as lRucinskil (l2006h did when 
establishing PL relation. W UMa which are known to have 
companions are marked in Table 1 . The 1 3 out 29 stars with 
o-jt/tt < 0.10 are marked, so after removing them, the 16 
W UMa systems left would definitely not satisfy statistical 
significance. We decided to include all systems up to the 
limit o-tt/tt < 0.15, which means uncertainties in the abso- 
lute magnitudes would be better than 0.33 mag. This process 
provided us with 3 1 W UMa stars after removing the ones 
with companions, which we think satisfy statistical signif- 
icance and accuracy for re-calibrating the PLC relation. W 
UMa stars with larger parallax errors are not included to 
avoid adding stars with less accurate absolute magnitudes. 



al 



20081) . which are 

Jo = J-0.887£;(B-y), (7) 

{J-H)o = {J-H)-0.322E{B-V), (8) 

iH-Ks)o^iH-Ks)-0.m3EiB-V). (9) 



In Table 2, LK corrected absolute magnitudes My and 
AIj for 3 1 W UMa stars were computed with most probable 
true parallax and adopted color excess E{B — V). The un- 
corrected absolute magnitudes [My {Hip) and Mj{Hip)] 
were also listed just for comparison. Fig. 3 compares cor- 
rected and uncorrected distances of the calibration sample. 
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Table 2 Calibration sample and data used in recalibration of the PLC relation for W UMa stars. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


Name 


7T 


CTtt/tT 


P 


SpT 


(B-V)o 


E(B-V) 


E(B-V) 


E(B-V) 


V™ax 


(B-V) 


E(B-V) 


J 


(J-H) 


(H-Ks) 


M(V) 


M(V) 


M(J) 


M(J) 




(mas) 




(day) 




INT^ 


OMI^ 


LIT-' 


SCH'' 


(mag) 


OBS= 


ADP'' 


(mag) 


(mag) 


(mag) 


(Hip) 


(true) 


(Hip) 


(true) 


eCiA 


33.32 


0.028 


0.59144070 


F2V 


0.35 


0.046 


0.000 


0.048 


4.74 


0.396 


0.046 ° 


4.052 


0.264 


-0.090 


2.218 


2.211 


1.632 


1.625 


V972 Her 


16.16 


0.038 


0.44309740 


F4V 


0.42 


0.046 


- 


0.123 


6.62 


0.466 


0.046 " 


5.785 


0.149 


0.103 


2.532 


2.519 


1.798 


1.786 


AEPhe 


20.36 


0.040 


0.36237274 


F8V 


0.53 


0.111 


0.000 


0.005 


7.56 


0.641 


0.000 ' 


6.577 


0.289 


0.107 


4.118 


4.104 


3.135 


3.121 


DOCha 


12.46 


0.045 


0.68144600 


F7V 


0.50 


0.032 


0.000 


0.024 


7.61 


0.535 


0.032 " 


6.634 


0.250 


0.067 


3.006 


2.988 


2.101 


2.083 


S Ant 


13.15 


0.053 


0.64834550 


A9V: 


0.30 


0.058 


0.010 


0.018 


6.40 


0.358 


0.010 ' 


6.087 


0.127 


0.095 


1.988 


1.963 


1.697 


1.673 


OUSer 


17.10 


0.055 


0.29676450 


F9/G0V 


0.56 


0.075 


0.010 


0.011 


8.10 


0.635 


0.010 '' 


6.946 


0.264 


0.112 


4.260 


4.234 


3.129 


3.102 


V759 Cen 


15.66 


0.059 


0.39399912 


F9V 


0.56 


0.034 


0.020 


0.017 


7.40 


0.594 


0.020 *> 


6.530 


0.242 


0.063 


3.342 


3.311 


2.517 


2.486 


GRVir 


18.53 


0.063 


0.34696950 


F7-8V 


0.50 


0.076 


0.010 


0.022 


7.80 


0.576 


0.010 *> 


7.048 


0.242 


0.044 


4.143 


4.108 


3.413 


3.378 


YYEri 


17.63 


0.067 


0.32150003 


G5V 


0.68 


-0.006 


0.000 


0.020 


8.10 


0.674 


0.000 '' 


7.031 


0.337 


0.106 


4.372 


4.331 


3.303 


3.262 


ISCMa 


9.80 


0.071 


0.61698000 


F3V 


0.38 


-0.008 


0.010 


0.012 


6.96 


0.372 


0.010 '' 


6.619 


0.151 


0.071 


1.931 


1.886 


1.612 


1.567 


CPHyi 


8.20 


0.072 


0.47940600 


FOV 


0.32 


0.027 


0.016 


0.019 


7.80 


0.347 


0.016"= 


7.187 


0.105 


0.099 


2.367 


2.320 


1.789 


1.743 


YYCrB 


11.10 


0.075 


0.37656400 


F8V 


0.53 


0.092 


0.010 


0.006 


8.48 


0.622 


0.010 '' 


7.640 


0.234 


0.097 


3.726 


3.675 


2.908 


2.857 


V757 Cen 


13.83 


0.078 


0.34316916 


F9V 


0.56 


0.084 


0.020 


0.022 


8.30 


0.644 


0.020 ' 


7.336 


0.304 


0.109 


3.996 


3.941 


3.077 


3.022 


FXEii 


13.33 


0.078 


0.29234500 


G9: 


0.78 


-0.007 


- 


0.005 


9.56 


0.773 


0.000 ° 


8.152 


0.407 


0.108 


5.239 


5.184 


3.831 


3.776 


RRCen 


9.46 


0.087 


0.60569200 


FOV 


0.32 


0.086 


0.030 


0.086 


7.27 


0.406 


0.086 " 


6.765 


0.113 


0.078 


1.951 


1.882 


1.636 


1.567 


GMDra 


9.84 


0.087 


0.33874120 


F5V 


0.45 


0.059 


0.030 


0.013 


8.66 


0.509 


0.030 '' 


7.736 


0.252 


0.043 


3.601 


3.533 


2.744 


2.675 


V535 Ara 


8.49 


0.101 


0.62930107 


A8V: 


0.27 


0.065 


0.020 


0.020 


7.17 


0.335 


0.020 '' 


6.989 


0.149 


0.021 


1.848 


1.753 


1.711 


1.616 


TYMen 


5.65 


0.106 


0.46166680 


A3/5V 


0.08 


0.174 


0.030 


0.087 


8.08 


0.254 


0.104" 


7.597 


0.088 


0.075 


1.623 


1.518 


1.370 


1.265 


V918 Her 


8.28 


0.107 


0.57481000 


A7V 


0.20 


0.077 


- 


0.128 


7.31 


0.277 


0.077 " 


6.703 


0.095 


0.013 


1,769 


1.662 


1.332 


1.225 


SXCrv 


10.33 


0.111 


0.31662090 


F6V 


0.48 


0.080 


0.040 


0.022 


8.99 


0.560 


0.040 ' 


7.927 


0.267 


0.075 


4.053 


3.937 


3.079 


2.963 


WYHor 


8.34 


0.113 


0.39894000 


G2IVA' 


0.63 


0.103 


- 


0.006 


9.39 


0.733 


0.103 " 


8.335 


0.321 


0.105 


3.798 


3.678 


2.971 


2.850 


EX Leo 


9.31 


0.113 


0.40860250 


F6V 


0.48 


0.049 


0.000 


0.015 


8.13 


0.529 


0.000 '' 


7.332 


0.217 


0.068 


3.095 


2.974 


2.297 


2.176 


OQVel 


5.07 


0.115 


0.58133800 


A3IV 


0.08 


0.100 


- 


0.116 


7.67 


0.180 


0.100° 


7.323 


0.037 


0.068 


1.010 


0.885 


0.884 


0.759 


NNVir 


8.90 


0.120 


0.48071484 


FO/IV 


0.32 


0.087 


0.000 


0.017 


7.60 


0.407 


0.000 "= 


7.060 


0.126 


0.114 


2.484 


2.347 


1.944 


1.807 


V2377 Oph 


9.46 


0.121 


0.42540100 


GO/ IV 


0.60 


0.056 


0.029 


0.022 


8.45 


0.656 


0.029 "= 


7.350 


0.275 


0.058 


3.380 


3.239 


2.344 


2.204 


VZPsc 


15.68 


0.123 


0.26125918 


K5V 


1.15 


0.120 


0.000 


0.026 


10.15 


1.270 


0.000 ■* 


8.158 


0.576 


0.176 


6.273 


6.127 


4.281 


4.135 


V351 Peg 


6.85 


0.125 


0.59329700 


A8V 


0.27 


0.074 


- 


0.025 


7.92 


0.344 


0.074 ° 


7.266 


0.099 


0.038 


2.019 


1.869 


1.529 


1.378 


TV Pic 


4.70 


0.125 


0.85198700 


A2V 


0.05 


0.092 


- 


0.010 


7.37 


0.142 


0.092 "^ 


7.096 


0.027 


0.072 


0.597 


0.446 


0.527 


0.376 


V870 Ara 


9.23 


0.134 


0.39978000 


F8 


0.53 


0.186 


- 


0.026 


8.80 


0.716 


0.026 = 


7.741 


0.264 


0.089 


3.722 


3.545 


2.720 


2.544 


V839 Cen 


10.93 


0.139 


0.33093400 


G2: 


0.63 


0.019 


- 


0.016 


8.82 


0.649 


0.019" 


8.315 


0.268 


0.146 


4.146 


3.955 


3.683 


3.492 


V386 Pav 


6.52 


0.141 


0.55184100 


A9V 


0.30 


0.010 


- 


0.022 


8.28 


0.310 


0.010° 


7.759 


0.101 


0.091 


2.517 


2.320 


2.019 


1.821 


1 : intrinsic cc 


lor, 2: Observed minus intrinsic colors, 3 
ss of observed minus intrinsic colors, b:[| 


talien from literature, 4: based on Schlegel et al. (1998), 5: Observed col 


3r, 6: Adopted color excess 
l.e:ISchle2eletalJ (19981) 






a: Color exce 


ucinslflH2006l), c:INordstrom et alj (20041), d:IRucinski 
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Fig. 3 Comparing corrected and uncorrected (Hipparcos) 
distances of W UMa stars with parallax errors (7-,^/'^ < 
0.15. 



3 Calibrations and discussions 

3.1 PLC relation according to uncorrected Hipparcos 
parallaxes 



There were 40 contact binaries used by lRucinski & Duerbeck 
(Il997h when establishing former PLC relation. Only 11 of 
them (e CrA, AE Phe, V759 Cen, V757 Cen, RR Cen, GR 
Vir, YY Eri, V535 Ara, TY Men, SX Crv, VZ Psc) are re- 
tained in our final calibration sample (Table 2). The others 
are eliminated either because they are members of multi- 
ple systems or their parallax errors are out of the acceptable 
limits of this study. 

In the first s tep, we were curious if t he coe fficients of the 
PLC relation by lRucinski & Duerbeckl(ll997h would change. 
Therefore, the uncorrected My of the W UMa stars in the 
calibration list are used in a regression analysis, which uses 
least squares fit for determining the three coefficients ap- 
pearing in the PLC relation. The regression analysis pro- 
duced similar coefficients. Therefore, the PLC relation of 



Rucinski & Duerbeckl (Il997h were kept as the reference re- 



lation with LK bias (uncorrected). 

Fig. 4a displays the absolute vis ual magnitudes which 



are pre dicted by the PLC relation of iRucinski & Duerbeck 



(Il997h against the absolute visual magnitudes from the Hip- 
parcos parallaxes directly. According to the regression anal- 
ysis, the correlation coefficient R^ = 0.916 and the stan- 
dard deviation from the diagonal is s = 0.35. 
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Fig. 4 a) Absolute visual magnitudes predicted by the 
PLC relation of Rucinski & Duerbeck (1997) versus abso- 
lute visual magnitudes of Hipparcos parallaxes, b) absolute 
visual magnitudes predicted by PLC relation corrected for 
LK bias versus absolute visual magnitudes corrected for LK 
bias. 



So, it can be concluded that having only more accurate 
parallax es (cr,r /tt < 0.15) than the orig inal sample {(Jt^I'k < 
0.24) of iRucinski & Duerbeck (Il997l) . the result does not 
justify re-calibrating the PLC relation, but only confirms an 
existing relation. 



the observed colors and adopted colorexcesses {E{B — V)), 
of the stars in Table 2, the coefficients of the linear equation 
M{V)NEW^alogP + b{B-V)o + c, (10) 

have been determined as a = -3.88(±0.62), b = 3.36(±0.34) 
and c — 0.0044(±0.14) by a regression analysis with a 
correlation coefficient R^ — 0.959 and standard deviation 
s — 0.27. According to the regression analysis it can be 
concluded that the internal error of the new relation is s.e. = 
s/Vn = 0.27/ \/3T = 0.048 mag. That is, in a case if peri- 
ods and intrinsic colors are errorless, then the standard error 
of predicted My is expected to be =F 0.048 mag. The stan- 
dard error of My should become larger if errors are intro- 
duced through orbital periods and intrinsic colors appearing 
in the relation. 

3.3 Comparing PLC relations with and without LK 
bias 

Comparing PLC relations of same format with and without 
LK bias may help to analyze and study the effect of LK bias 
on the PLC relation. There could be two basic types of com- 
parison: The two relations may be compared on their own 
grounds to see how good their predicting powers are ac- 
cording to their own conditions; or they could be compared 
to see what the difference is in predicting true absolute mag- 
nitudes. 

Consequently, for the first type of comparison, the 
M{V)new points are plotted against the corrected abso- 
lute magnitudes coming from the corrected parallaxes in 
Fig. 4b. Although, the correlation coefficient is shghtly bet- 
ter (i?2 = 0.928) and the standard deviation s = 0.32 is 
smaller, the apparency of Fig. 4b is not much different Fig. 
4a. Both figures are consistent with each other. Thus, LK 
correction is justified. 

The difference between the two PLC relations can be 
observed better in Fig. 5 where the predictions of both rela- 
tions are plotted against the corrected absolute magnitudes. 
Although both data sets remain within the 2a reliability 
limit, it is clearly shown in F ig. 5 that the PLC relation of 
Rucinski & D uerbeckl (ll997h systematically underestimates 
the absolute brightness of the same sample stars. 

A similar comparison can be done on a distance scale 
rather than a magnitude scale. Fig. 6 displays su ch a com- 



pariso n. It is clearly shown that the PLC relation of lRucinski & Duerbeck 
(Il997h underestimates the distances with respect to the pre- 
dictions of the PLC relation from this study especially when 
the distances are large. If the distances are smaller than 100 
pc, apparently, the difference becomes negligible. This must 
be because; LK corrections become noticeable if (Jt^/tt > 
0.05 and increases as (Jt^/tt increases. 



3.2 PLC relation according to LK bias corrected 
parallaxes 

Using the corrected absolute magnitudes (My), orbital pe- 
riods (P) and the intrinsic {B — V)o colors determined from 



3.4 Calibrating PLC relation using 2MASS data 

PLC relation of W UMa stars could be calibrated using dif- 
ferent color indices and color excesses. A PLC relation us- 



ing {V — I)o already exists (IRucin ski. 2004). The advan- 
tage of shifting PLC relation towards the infrared colors is 
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Fig. 5 Absolute visual magnitudes predicted by two PLC 
relations (Rucinski & Duerbeck 1997 and this work) plot- 
ted against the absolute visual magnitudes directly from 
Pogson's relation using LK bias corrected parallaxes. The 
dashed-dotted and dashed lines represent prediction limits 
for ±1(7 and ±2cr, respectively. 
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Fig. 6 The distances according to both PLC relations are 
plotted against the LK bias corrected distances. 



obvious; extinction is less, so uncertainty of interstellar ab- 
sorption becomes relatively reduced. 

The method of calibration is the same. Absolute mag- 
nitudes, and then intrinsic colors through (8) and (9) were 
computed using the data in Table 2 for the same sample (3 1 
stars). Finally a regression analysis is done for determining 
the numerical coefficients. It has been tried for many dif- 
ferent forms of the relation such as; with one color or both 
(J — H)q and {H — Ks)o and even to include {B — V)q 
color besides the terms of log P including a free constant. 
The highest correlation have been found with the following 

M{J) = alogP + b{J - H)o + c{H - Ks)q + d, (11) 

where a = -3.17(±0.62), h = 3.84(±0.62), c = L81(±1.16) 
and d = 0.24(±0.14), which gives a correlation coefficient 
E? — 0.927 and a standard deviation s = 0.26 mag. Con- 
sequently, the internal error of the relation is ^0.047 mag. 
Since best correlation was achieved by using two rather than 
a single color term, it seems using two color terms increases 
the correlation. This is expected because each color term 
contributes to improving the correlation. The constant term 
which has an absolute value comparable to its uncertainty 
may not be real so it can be avoided in the relation. 

The predicted M{J) values by Eq. 1 1 are plotted against 
LK corrected absolute magnitudes (Fig. 7). Despite, 2MASS 
data are single epoch observations, which were not guaran- 
teed to be at the maximum brightness phase of the systems, 
the calibration gave reliable relation which can be used in 
predicting the true distances of W UMa stars from their pe- 
riods and infrared colors. It is even more interesting that the 
standard deviation of data and the internal error of the PLC 
relation in infrared are slightly smaller than the standard de- 
viation and internal error of the PLC relation in the visual 
although the PLC relation is slightly less correlated at in- 
frared. 

The distance predictions of the PLC relation are com- 
pared to LK bias corrected distances in Fig. 8. The internal 
error of the relation is =f0.022. The PLC relation at infrared 
could have been improved more if maximum brightnesses 
at these infrared colors were available. 



4 Conclusions 

In this study, a statistical bias, which is classically known 
as Lutz-Kelker (LK) bias, has been introduced to the PLC 
relation of W UMa stars. Since parallax errors remain after 
a LK bias correction, studies such as ours, should still be 
careful in setting up an upper limit on the reliability of input 
parallax data. In this study, we have preferred to work with 
parallax data more accurate than 15%. 

Using LK corrected parallaxes (or distan ces), we have 



re-calibrated the classical PLC relation of Rucinski & Duerbeck 



(Il997h . The new relation has been compared with the classi- 
cal relation. It has been found that old and new relations are 
nearly equivalent in predicting a distance up to about 100 
parsecs. The difference between the predictions increases 
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Fig. 7 Predicted M(J) are compared to LK corrected ab- 
solute magnitude. The dashed-dotted and dashed lines rep- 
resent prediction limits for ilcr and ±2cr, respectively. 
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Fig. 8 Distances predicted by the PLC relation in the in- 
frared are compared to the LK corrected distances. 




and may become intolerable as the distances of stars in- 
crease. The new relation is valid in the ranges (0.26 < P{day) < 
0.87, 0.04 <{B- V)o < L28 and 0.43 < M(y) < 6.55). 
We have also produced a PLC relation which can be 
used by 2MASS data. Despite 2MASS data are single epoch 
observations, which are not guaranteed to be at the max- 
imum brightness phase of the W UMa stars in the cali- 
bration sample (Table 2), the established relation has been 
found surprisingly consistent and reliable to predict LK cor- 
rected distances. This new relation is valid in the ranges 
(0.26 < Piday) < 0.87, -0.01 < (J - H)o < 0.58, 
-0.10 < {H - Ks)o < 0.18 and 0.54 < M{J) < 4.63). 
Observers are encouraged to obtain light and color curves of 
W UMa stars at 2MASS colors. If the apparent magnitudes 
at maximum brightness phases of W UMa stars are provided 
with 2MASS colors, the current PLC relation using infrared 
colors could be improved to provide more reliable distances. 
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